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Abstract 

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Psyllidae), is an introduced 
pest in Mexico and a vector of huanglongbing, a lethal citrus disease. Estimations of the habitat 
distribution and population growth rates of D. citri are required to establish regional and area- 
wide management strategies and can be used as a pest risk analysis tools. In this study, the habitat 
distribution of D. citri in Mexico was computed with MaxEnt, an inductive, machine-learning 
program that uses bioclimatic layers and point location data. Geographic distributions of devel- 
opment and population growth rates were determined by fitting a temperature-dependent, non- 
linear model and projecting the rates over the target area, using the annual mean temperature as 
the predictor variable. The results showed that the most suitable regions for habitat of D. citri 
comprise the Gulf of Mexico states, Yucatan Peninsula, and areas scattered throughout the Pacif- 
ic coastal states. Less suitable areas occurred in northern and central states. The most important 
predictor variables were related to temperature. Development and growth rates had a distribution 
wider than habitat, reaching some of the northern states of Mexico. Habitat, development, and 
population growth rates were correlated to each other and with the citrus producing area. These 
relationships indicated that citrus producing states are within the most suitable regions for the oc- 
currence, development, and population growth of D. citri, therefore increasing the risk of 
huanglongbing dispersion. 

Resumen 

El psilido asiatico de los citricos, Diaphorina citri Kuwayama (Hemiptera: Psyllidae) es una 
plaga introducida en Mexico, y es vector del huanglongbing, una enfermedad letal de los citricos. 
La estimacion de la distribucion del habitat y tasas de crecimiento poblacional de D. citri se 
requieren para establecer una estrategia de manejo a nivel regional y como herramienta de 
analisis de riesgo. En este estudio, la distribucion del habitat de D. citri en Mexico se calculo con 
MaxEnt, un programa inductivo, de aprendizaje-maquina que emplea capas bioclimaticas y datos 
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puntuales de localizacion. Tambien se determinaron las distribuciones geograficas de las tasas de 
desarrollo y crecimiento poblacional por el ajuste de un modelo no-lineal, dependiente de la 
temperatura y proyectando las tasas sobre el area objetivo usando la temperatura media anual 
como variable predictora. Los resultados mostraron que las regiones mas apropiadas para el 
habitat de D. citri comprenden los estados del Golfo de Mexico, Peninsula de Yucatan y areas 
dispersas de los estados costeros del Pacifico. Areas menos apropiadas ocurren en los estados del 
centro y norte. Las variables predictoras mas importantes estuvieron relacionadas con la 
temperatura. Las tasas de desarrollo y crecimiento poblacional presentaron una distribucion mas 
amplia que el habitat, alcanzando algunos de los estados del norte. El habitat y las tasas de 
desarrollo y crecimiento poblacional estuvieron correlacionadas entre si y con la superficie 
plantada con citricos. Estas relaciones sugieren que los estados productores de citricos se 
encuentran dentro de las areas mas apropiadas para la ocurrencia, desarrollo y crecimiento 
poblacional de D. citri, por tanto incrementando el riesgo de dispersion del huanglongbing. 
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Introduction 



The Asian citrus psyllid, Diaphorina citri 
Kuwayama (Hemiptera: Psyllidae) is an im- 
portant pest of commercial citrus worldwide. 
It has a variable life cycle, and egg to adult 
development is temperature dependent (Liu 
and Tsai 2000), as is the population growth 
rate, as measured by the intrinsic rate of in- 
crease (rm) (Nava et al. 2010). Diaphorina 
citri feeds on shoots, leaves, and buds of cit- 
rus plants (Yang et al. 2006). Also, it vectors 
the pathogen associated with huanglongbing, 
a lethal citrus disease (da Graca and Korsten 
2004; Bove, 2006). The pest has a worldwide 
distribution across tropical and subtropical 
countries (Yang et al. 2006). In Mexico, D. 
citri was detected in 2002 (Lopez-Arroyo et 
al. 2009) and it has dispersed throughout the 
country. Its presence has been reported in all 



of the citrus producing states (Lopez-Arroyo 
et al. 2009). 

Huanglongbing is a bacterial disease caused 
by the phloem-limited a-proteobacteria Can- 
didates Liberibacter spp (Rhizobiales: 
Rhizobiaceae). Three species of Candidatus 
have been determined: Ca. L. asiaticus, Ca. L. 
africanus, and Ca. L. americanus (Bove 
2006). Huanglongbing threatens the citrus in- 
dustry because it causes the death of the trees 
and is difficult to control. The disease is con- 
sidered a serious factor affecting the 
production of commercial citrus orchards, it 
has no cure, and only preventive measures can 
deter the spread of the disease once it has been 
established. Huanglongbing occurs in coun- 
tries in Africa, Asia, and America (Bassanezi 
et al. 2010). It is present in Brazil, U.S.A., 
Central America, Caribbean Islands, and Mex- 
ico (National Research Council 2010). In 
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Brazil, huanglongbing caused the removal of 
thousands of trees, and producers have allo- 
cated economic resources to reduce damage 
and consequently decrease the dispersion 
(Bassanezi et al. 2010). In Mexico, the disease 
was first detected in Yucatan in 2009, and af- 
ter that it appeared in some municipalities 
across the Pacific coastal states (SENASICA 
201 1). Regulatory control measures have been 
undertaken to control the disease, and a na- 
tionwide alert campaign has been launched, 
mainly in citrus producing states (Mayorga 
2010). Mexico ranks among the five top citrus 
producing countries worldwide (FAO 2008), 
so the occurrence of D. citri and huanglong- 
bing is a serious threat to the Mexican citrus 
industry (Salcedo et al. 2010). 

The geographic distribution of introduced 
pests is an important component for risk anal- 
ysis and making decisions on control methods. 
Species distribution modeling is a key tool to 
study the capacity of a species to occur in a 
given region (Elith and Leathwick 2009). Es- 
timates are based on deductive, inductive, and 
hybrid methods. Deductive methods compute 
species distribution based on climatic re- 
quirements and do not require presence data. 
Inductive methods rely on actual sampling 
points, thus providing an estimate based on 
presence of the pest. Hybrid methods combine 
both (Venette et al. 2010). MaxEnt is an in- 
ductive, machine-learning program that 
computes the habitat or niche distribution for 
a given organism (Phillips et al. 2006) and has 
performed well when compared to other 
methods (Ortega-Huerta and Peterson 2008). 
However, for a large sample size, boosted de- 
cision trees performed slightly better than 
MaxEnt, but they did not differ significantly 
from each other (Wisz et al. 2008). Like other 
distribution models, MaxEnt uses environ- 
mental layers and presence-only data to 
predict habitat distribution for a given organ- 
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ism (Phillips et al. 2006); point occurrence 
data are most commonly generated, as com- 
pared to presence-absence data (Tsoar et al. 
2007; Elith et al. 2011). 

The habitat distribution gives a view of where 
a species could establish, but development and 
population growth rates are biophysical pro- 
cesses that provide additional information on 
how insect populations respond to abiotic fac- 
tors (Buckley et al. 2010). Most studies have 
focused on studying the effect of temperature 
on these rates, because insects are poikilo- 
thermic organisms (Wagner et al. 1984). The 
development and population growth of D. citri 
has been measured on different citrus species, 
including the orange jasmine, Murray a panic- 
ulata, and degree-days models have been 
computed to describe development (Nava et 
al. 2010). However, the development of in- 
sects is intrinsically non-linear (Wagner et al. 
1984), and, to account for the variability in 
temperatures occurring at country level, this 
non-linearity should be included. Several non- 
linear models have been proposed and used to 
represent temperature-dependent variables, 
which range from three- to six-parameter 
models (Wagner et al. 1984). Briere 1 model 
has been compared to more complex models 
with similar performance; nonetheless, it uses 
only three parameters, allowing it to handle 
small datasets (Shi and Ge 2010), and it has 
been applied to model aphid development, 
among other organisms (Jalali et al. 2010). In 
addition to habitat prediction, an estimate of 
the distribution of potential growth and devel- 
opment may shed light on the biological 
capacity of D. citri to thrive in the different 
agroclimatic regions in Mexico. Given that 
huanglongbing is already present in Mexico, 
the objective of this study was to know the 
potential distribution of habitat, development, 
and population growth rates of D. citri in this 
country. 
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Materials and Methods 

Habitat distribution 

Habitat distribution of D. citri was estimated 
with MaxEnt version 3.3.3 (Phillips et al. 
2006). MaxEnt generates a variable that 
measure the probability to occur in a given 
region and is considered as the realized niche 
(Phillips et al. 2006); probability ranges from 
0 (least likely) to 1 (most likely) to occur. The 
model requires two types of data: species 
point location data and bioclimatic layers. Di- 
aphorina citri point location data were 
collected and geo-referenced by personnel of 
the Mexican Plant Health State Committee 
during 2008 and 2009 across most of the cit- 
rus producing states as a result of a 
nationwide huanglongbing detection effort 
(Robles-Garcia and Delgadillo-Villanueva 
2008). Bioclimatic data were retrieved from 
Atlas Digital de Mexico (Fernandez-Eguiarte 
et al. 2010), comprising 19 raster layers de- 
rived from temperature and precipitation and 
representing extreme values, seasonality, and 
annual trends (Table 1). These layers were 
curated from the WorldClim database 
( www.worldclim.org ), were derived from his- 
toric global source data ranging from 1950 to 
2000, and had a spatial resolution of 30" or 
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Table 1. Bioclimatic variables used to fit the Maxent model. 


CODE 


Bioclimatic Variable 


BIO-1 


Annual Mean Temperature 


BIO-2 


Mean Diurnal Range (Mean of monthly (max temp - min temp)) 


BIO-3 


Isothermality (BIO-2/BIO-7) (* 100) 


BIO-4 


Temperature Seasonality (standard deviation * 100) 


BIO-5 


Max Temperature of Warmest Month 


BIO-6 


Min Temperature of Coldest Month 


BIO-7 


Temperature Annual Range (B10-5-BIO-6) 


BIO-8 


Mean Temperature of Wettest Quarter 


BIO-9 


Mean Temperature of Driest Quarter 


BIO- 10 


Mean Temperature of Warmest Quarter 


BIO- 11 


Mean Temperature of Coldest Quarter 


BIO- 12 


Annual Precipitation 


BIO- 13 


Precipitation ofWettest Month 


BIO- 14 


Precipitation of Driest Month 


BIO- 15 


Precipitation Seasonality (Coefficient of Variation) 


BIO- 1 6 


Precipitation ofWettest Quarter 


BIO- 17 


Precipitation of Driest Quarter 


BIO- 18 


Precipitation of Warmest Quarter 


BIO- 19 


Precipitation of Coldest Quarter 
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about 0.86 km 2 at the Equator; these layers 
were generated by interpolation with a thin- 
plate smoothing spline algorithm (Hijmans et 
al. 2005). The fitness of the model was meas- 
ured by reserving 30% of sampling points for 
testing purposes (test data) and computing the 
area under the curve (AUC) to compare the 
predictive power of the model when using the 
test data against the training data (Phillips and 
Dudik 2008). AUC values less than 0.7 indi- 
cate a low performance model, values 
between 0.7 and 0.9 indicate a useful model, 
and values higher than 0.9 indicate a very 
good model (Baldwin 2009). The importance 
of the bioclimatic variables was estimated 
with jackknife where the relative contribution 
of each variable was evaluated either as its 
single contribution or by removing it from the 
predictor variable set (Phillips et al. 2006). 
For the rest of the model parameters, the de- 
fault values were used, as recommended by 
Phillips and Dudik (2008). The citrus growing 
area was not used as a bioclimatic layer due to 
its low resolution; instead, a regression analy- 
sis at state level was performed and described 
later. 

Development and population growth rate 
distributions 

The geographical potential for development 
rate, which is the reciprocal of development 
time (1/d) and population growth rate (rm), 
was computed in a different way as for habitat 
and applied a biophysical model that describes 
these indices as a function of temperature. A 
two-step approach was used: first, the parame- 
ters of the Briere 1 model (Briere et al. 1999) 
were estimated as a function of temperature 
using the NLIN procedure with SAS v. 9.1 
(SAS Institute 2004). The model is: 



f(T) = a(T-x mm )^ max -T) 
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where f(T) is either development rate (1/d) or 
the intrinsic rate of increase (rm) as a function 
of temperature (T). The model parameters are: 
a, an empirical constant and x m in and x max , 
which represent the minimum and maximum 
threshold temperatures; T is the constant tem- 
perature at which individuals are kept to 
measure their development. 

Development rate data, from egg to adult, 
were used to fit the model and were obtained 
from Nava et al. (2010). For the intrinsic rate 
of increase, data published by Liu and Tsai 
(2000) were analyzed. In the second step, 
once the model parameters were estimated (a, 
Xmin, and Xmax), the geographic distribution of 
1/d and rm was calculated by using the aver- 
age annual temperature raster layer for the 
whole country as the predictor variable T 
(Fernandez-Eguiarte et al. 2010) and was used 
as input to the Briere 1 model to obtain the 
corresponding response f(T), that is, T is the 
raster layer containing the average annual 
temperature and a, x m i n , and x max substituted 
by their nonlinear estimates. These map alge- 
bra operations were done with Octave version 
3.2 (Eaton et al. 2008). Finally, to determine 
the similarity among the three suitability indi- 
ces (probability, 1/d, and rm) and its 
relationship with the citrus growing area (ha) 
for year 2009 (SIAP 2010), three linear mod- 
els of the form Y = Po + PiX were estimated at 
state level; therefore, average distribution val- 
ues were computed at this level. For each 
distribution, the state area was clipped with 
the corresponding state polygon, and then the 
average value was computed. These opera- 
tions were performed with SAGA version 
2.0.7 (SAGA 2011). For the relationship be- 
tween probability and the citrus growing area, 
mean values were transformed to natural loga- 
rithms (In). The model parameters f3o and Pi 
were estimated by regression analysis using 
PROC REG in SAS version 9.1 (SAS Institute 
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2004). All the maps were composed by over- 
laying the distribution raster layers with the 
political boundaries with ArcGis version 9.3.1 
using a WGS84 projection (ESRI 2009); the 
state-level vector layer was obtained from 
INEGI (INEGI 2011). 

Results 

Habitat distribution of Diaphorina citri 

The presence of D. citri in Mexico was de- 
tected for most of the citrus producing states 
(Figure 1) and reflected the effort to sample 
the most important citrus growing areas, such 
as Colima, Michoacan, and Veracruz, and the 
potential huanglongbing entry points such as 
the northern border states (Baja California, 
Nuevo Leon, Tamaulipas) and the Yucatan 
peninsula. The estimated probability of habitat 
distribution of D. citri is shown in Figure 2. 
High values of presence occurred in parts of 
Baja California Sur, Sonora, Sinaloa, Colima, 
Michoacan, and Guerrero, along with the Gulf 
of Mexico states and the Yucatan peninsula. 
Furthermore, northern and central states like 
Chihuahua, Coahuila, Durango, Zacatecas, 
Aguascalientes, and Guanajuato had the low- 
est values. In general, D. citri had a high, 
uniform distribution in all of the Gulf of 
Mexico states and was scattered throughout 
the Pacific states. 

The model fitness for the estimated D. citri 
distribution was measured with the AUC 
(Figure 3). The red line indicates the fit with 
the training data, and the blue line indicates 
the actual model performance with the inde- 
pendent test data. Both were similar in their 
predictive value, and model prediction is sta- 
tistically different from a random one (p < 
0.01, binomial test). AUC values over 0.7 are 
considered appropriate for predictive purposes 
(Phillips and Dudik 2008; Baldwin 2009); 
therefore, the model was deemed appropriate 
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to show where D. citri can occur. MaxEnt us- 
es several bioclimatic layers for predictive 
purposes; however, their relative importance 
differs among them. The contribution of the 
predictor variables to model fitting (train data) 
is shown in Figure 4. The five most important 
by their individual contribution (blue bars) 
were temperature-related, two refer to extreme 
cold (BIO-6, BIO-11), and the others to aver- 
age (BIO-1) or fluctuating values (BIO-4, 
BIO-7); the less important variables were re- 
lated to high temperatures and to precipitation. 
The same figure shows that when a variable is 
excluded from the model (green bars), it con- 
tains less valuable information than by itself, 
that is, the whole model is more resilient to 
variable exclusion. On the other side, Figure 5 
shows variable contribution for the test data 
and has similar values as those for model fit- 
ting. Overall, the most important variables 
have critical information either to fit the mod- 
el (train data) or to predict values (test data). 

Development and population growth rate 
distributions 

The Briere 1 model was significant for devel- 
opment (F 3 , 3 = 447.79; p < 0.01). The 
estimated parameters were a = 4.5E-5, x m i n = 
12.6° C, and x max = 38.1° C. When this model 
was applied using the annual mean tempera- 
ture raster layer T as the input, the 
development rate had a distribution wider than 
niche, presenting high values along most of 
the coastal states, the Yucatan peninsula, and 
Baja California Sur (Figure 6). This result in- 
dicates that D. citri had a high potential for 
development in most of the suitable locations. 
MaxEnt estimates were constrained by the 
presence-only records, while the computations 
of (1/d) were constrained only by the lower 
and upper threshold temperatures, thus the 
wider distribution. 
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Regarding the potential of population growth, 
the parameters for the Briere 1 model were a 
= 2.3E-4, „ min = 13.2° C and „ max = 31.4° C, 
and the effect of temperature was significant 
(F 3 ,2 = 61.83; p = 0.016). This model was 
solved over the annual mean temperature T, 
and the resultant geographic distribution of rm 
is shown in Figure 7. Like the distribution of 
development rate, the higher values occurred 
along the coastal states of the Pacific Ocean 
and the Gulf of Mexico, including the Yuca- 
tan peninsula and Baja California Sur. 
Therefore, these regions presented a high po- 
tential for the population to grow. 

At the state level, the three indices were relat- 
ed. First, development rate and habitat had a 
significant linear relationship (Fi j2 6 = 103.4,/? 
< 0.01; r 2 = 0.79) (Figure 8), hence the simi- 
larity in their distribution maps. Second, the 
relationship between population growth and 
development rates showed a closer similarity 
(Fi j2 6 = 2947.7,/? < 0.01; r 2 = 0.99) (Figure 9) 
because both are affected in a similar way by 
temperature. Third, the effect of citrus host on 
habitat was estimated by the relationship be- 
tween habitat distribution and the citrus 
growing area, and this was significant (Fi, 2 6 
31.3./; < 0.01) (Figure 10). In this case, how- 
ever, the trend was coarse, given the medium 
value of r 2 (0.54). This model indicated that 
the larger the growing area, the higher the 
probability of pest presence. 

Discussion 

Habitat distribution of Diaphorina citri 

Diaphorina citri was first detected in Mexico 
in 2002 and it has extended to most of the 
country (Lopez-Arroyo et al. 2009). It is 
thought to be native from tropical and sub- 
tropical regions of Asia (Liu and Tsai 2000; 
Halbert and Manjunath 2004; Yang et al. 
2006). The habitat distribution estimated with 
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MaxEnt roughly agreed with the biogeograph- 
ical region known as the neotropical zone 
(Morrone 2002), which ranges from Central 
America to central and southern Mexico, and 
was consistent with the distribution of other 
psyllid species, such as Heteropsylla cubana 
(Hodkinson 1989). However, some of these 
regions had low probability values, such as 
Guerrero, Oaxaca, and Chiapas. Host occur- 
rence is an important factor affecting D. citri 
distribution but was not included in the 
MaxEnt analysis; in fact, Guerrero, Oaxaca, 
and Chiapas have small areas planted with 
citrus (Figure 1), thus the chance to find D. 
citri was reduced. On the contrary, northern 
states like Nuevo Leon and Tamaulipas, 
which are beyond the neotropical region, had 
high probability values and large areas planted 
with citrus (SIAP 2010). Therefore, the distri- 
bution of D. citri depends on the availability 
and distribution of citrus plants as well (Ro- 
bles-Garcia and Delgadillo-Villanueva 2008). 
This factor was accounted for by the relation- 
ship between probability values and the 
surface planted with citrus, as modeled by the 
linear model (Figure 7), and was significant, 
as expected. 

Moreover, the main variables affecting proba- 
bility distribution were related to temperature, 
from which the minimum temperature of 
coldest month seemed to be the most im- 
portant; in consequence, it appears that 
presence was most affected by cold tempera- 
tures. These findings agree with those 
reported by Yang et al. (2006), who cite am- 
bient temperature as the most important factor 
limiting the distribution of D. citri, and partic- 
ularly average temperature of the coldest 
month, minimum temperature, and the dura- 
tion of cold weather. Extreme temperatures 
have long been considered to affect plant dis- 
tribution (Magarey et al. 2008), and, for 
insects, non-linear development models ex- 
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plicitly account for minimum and maximum 
temperatures that constrain development and 
therefore occurrence of a given insect species 
(Wagner et al. 1984). This was explicitly ac- 
counted for by the Briere 1 model, which 
restricts development below 12.6° C. Other 
estimates of lower temperature threshold for 
degree-day models indicated similar values, 
ranging from 11.1 to 14° C (Liu and Tsai 
2000; Nava et al. 2010). Oviposition of D. cit- 
ri is another process that extreme 
temperatures affect. It has been demonstrated 
that D. citri does not oviposit below 16° C or 
above 41.6° C (Hall et al. 2011). These previ- 
ous biological studies support the importance 
of considering extreme temperatures as factors 
influencing D. citri distribution. 

Regarding performance of MaxEnt, the com- 
puted AUC values for both test and train data 
showed a good fit, given they exceeded the 
acceptable threshold of 0.7 (Phillips and Du- 
dik 2008; Baldwin 2009). On the other hand, 
the sample size used in this work comprised 
4,514 records for training and 1,934 for test- 
ing and, given that MaxEnt was considered 
among the best with predictive power across 
different sample sizes (Hernandez et al. 2006; 
Wisz et al. 2008), we expected the estimated 
habitat distribution to be representative of the 
insect's potential to occur in Mexico. This 
sample size was larger than the range of 10 to 
30 point data, which is considered as too re- 
strictive to make inferences on the species 
distribution (Wisz et al. 2008), and exceeds 
the 50-100 points estimated to obtain a useful 
distribution model (Wisz et al. 2008; Franklin 
and Miller 2009). On the other hand, the best 
predictor variables are the same for testing 
and for training the model, therefore indicat- 
ing that they are both valuable to modeling as 
for predicting. In any case, the best predictor 
variables indicated that extreme temperatures 
affect D. citri distribution, which is supported 
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by independent research (Yang et al. 2006; 
Hall et al. 2011). 

Development and population growth rate 
distributions 

The distributions of development and popula- 
tion growth rates were wider than for MaxEnt 
estimates. In the case of rm distribution, it was 
very similar to development, which probably 
results from the fact that as the development 
rate increases, the generation time decreases 
and rm increases (Poole 1974). The wider dis- 
tribution of both variables as compared to 
habitat seems to occur because models based 
on underlying physiological or population re- 
sponses, like development rate and rm, 
usually estimate the larger, potential niche, as 
compared to the smaller, realized niche esti- 
mated by point data (Buckley et al. 2010; 
Venette et al. 2010). Another probable factor 
affecting a wider distribution is that these in- 
dices were computed using a single predictor 
variable, while habitat was constrained by ac- 
tual presence and by several bioclimatic 
layers. Both rates were dependent on annual 
mean temperature only, and additional envi- 
ronmental factors influence the inherent 
capacity to develop and grow; however, tem- 
perature is considered as the most important 
variable that affects insect development and 
has been considered as the primary factor to 
estimate species distribution when mechanis- 
tic models are used (Buckley et al. 2010). For 
example, in Australia, temperature was select- 
ed as the main factor affecting D. citri 
distribution under climate changes (Auram- 
boutetal. 2009). 

The close relationship among probability val- 
ues and citrus host area was most likely 
because citrus plants are the preferred host of 
D. citri (Tsai and Liu 2000; Halbert and 
Manjunath 2004). In general, the three linear 
models suggest that habitat was related to cit- 
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rus area, and states with the largest area plant- 
ed with citrus had a high potential for D. citri 
to establish and grow; therefore, the risk of 
huanglongbing transmission might increase. 
The presence of huanglongbing in Mexico 
occurs in regions that, according to these re- 
sults, can be considered as among the best 
suited for the vector to develop: Baja 
California Sur, Chiapas, Colima, Michoacan, 
Sinaloa, Hidalgo, Nayarit, Quintana Roo, San 
Luis Potosi, and Yucatan (SENASICA 2011). 
These results highlight the necessity to im- 
plement immediate measures aimed to detect, 
eliminate, confine, or control the huanglong- 
bing inoculum sources. Important citrus 
producing areas (Figure 1) that have not been 
invaded by huanglongbing are the states of 
Veracruz, Tamaulipas, and Nuevo Leon. In 
such places, the need to establish a monitor- 
ing-alarm system for the rapid detection and 
actions against the disease is crucial. 

Overall, the results indicated that coastal 
states and Yucatan and Baja California penin- 
sulas are highly suitable for the occurrence, 
fast development, and population growth of 
D. citri. Therefore, these regions can be con- 
sidered of high risk for huanglongbing 
transmission given the potential of the vector 
to establish and increase its populations. The- 
se results also suggest that eradication of D. 
citri would be difficult given the suitability for 
its development. As in other countries, 
measures to contain the spread of huanglong- 
bing should be effected by the integration of 
different methods, for example, area-wide 
management of the vector (Qureshi and 
Stansly 2010), use of healthy plant material, 
and removal of infected plants (Bassanezi et 
al. 2010; Mayorga 2010; National Research 
Council 2010). 
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Figure 4. Jackknife test of variable importance for train data: individual variable contribution (blue bar), contribution when a 
given variable is excluded (green bar), whole set of variables (red bar). Variable names are listed in Table I. High quality figures 
are available online. 
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Figure 5. Jackknife test of variable importance for test data: individual variable contribution (blue bar), contribution when a 
given variable is excluded (green bar), whole set of variables (red bar). Variable names are listed in Table I . High quality figures 
are available online. 
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